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In this second paper, we develop translerable semi-empirical parameters for the technologically 
important material, silicon, using Extended Hiickel Theory (EHT) to calculate its electronic struc- 
ture. The EHT-parameters are optimized to experimental target values of the band dispersion of 
bulk-silicon. We obtain a very good quantitative match to the bandstructure characteristics such as 
bandedges and effective masses, which are competitive with the values obtained within an sp 3 d 5 s* 
orthogonal-tight binding model for silicon—. The transferability of the parameters is investigated 
applying them to different physical and chemical environments by calculating the bandstructure 
of two reconstructed surfaces with different orientations: Si(100) (2x1) and Si(lll) (2x1). The 
reproduced n- and 7r*-surface bands agree in part quantitatively with DFT-GW calculations and 
PES/IPES experiments demonstrating their robustness to environmental changes. We further ap- 
ply the silicon-parameters to describe the ID band dispersion of a unrelaxed rectangular silicon 
nanowire (SiNW) and demonstrate the EHT-approach of surface passivation using hydrogen. Our 
EHT-parameters thus provide a quantitative model of bulk-silicon and silicon-based materials such 
as contacts and surfaces, which are essential ingredients towards a quantitative quantum transport 
simulation through silicon-based heterostructures. 



I. INTRODUCTION 

Silicon is the dominant component in fabrication of 
semi-conducting devices and continues to play a key 
role in future nanoelectronic devices. Recent STM- 
experiments of molecules to highly doped silicon (n- or 
p-typc doping) showed NDR-behaviour at room temper- 
ature in the molecular current-voltage characteristics^ 
which might be a modest initial step towards a molecule 
based electronics. The NDR-behavior, however, was pre- 
dicted and theoretically demonstrated by Datta and co- 
workers using a physical resonable model for the elec- 
tronic structure of silicon. Therefore, an correct model 
for the electronic structure of all constituents is essen- 
tial for a quantitative modeling of quantum transport on 
nano- and atomistic scales. 

DFT-based approaches are well-suited to determine 
electronic and atomic material properties which depend 
on the total energy^ Abinitio methods in various approx- 
imations such as DFT-LDA/GGA have been successful 
to describe properties of molecules and metalsj^ but 
they are less benchmarked for the electronic properties 
of semi-conducting materials, particularly for multiple 
bandedges and effective masses - both of them are impor- 
tant components for a quantitative simulation of quan- 
tum transport. For semi-conductors one well known fail- 
ure of DFT-LDA/GGA calculations is that the bandgap 
is systematically underestimated by about a factor of 2^ 
which makes a quantitative modeling for example of sil- 
icon difficult. Quantitative correct bandgaps can be ob- 
tained within the GW-approximation,® which is compu- 
tationally expensive and might thus have limited usage 
at least for transport through large nanostructures. 

At the other end arc less rigorous, but computationally 



less expensive semi-empirical methods such as orthogonal 
and non-orthogonal tight binding approaches. Here, the 
electronic structure is not rigorously calculated, but de- 
termined by an optimization of the free parameters such 
as the matrix elements of the Hamiltonian and overlap 
to match bandedges and effective masses. Orthogonal 
tight-binding (OTB) approaches have been extensively 
developed in the past and further optimized to describe 
the electronic structure of bulk semi-conductors such as 
silicon and germanium, for exampl o 7 i 8 i 9 . However, the 
transferability of the bulk optimized parameters - a prob- 
lem of any scmi-cmipirical approach - has not been clearly 
demonstrated, for example, by calculating the electronic 
structure of (re-constructed) surfaces of different orien- 
tations. In turn, their robustness has been mainly shown 
for atomic structure optimization of bulk, surfaces, and 
finite size clusters, which requires the calculation of the 
total energy. 10 i 11 i 12 i 13 i 14 i 15 i 16 The latter, however, is just 
an integral property of the band dispersion and eventual 
errors in the bandstructure cancel^ 

In this paper we use a non-orthogonal tight-binding 
scheme - Extended Hiickel Theory (EHT) - to calcu- 
late the electronic structure of the technologically im- 
portant material silicon. In our first paper, we applied 
the EHT approach to carbon nanotubes demonstrating 
the transferability of the EHT-parameters for carbon to 
small diameter tubes as well as to a strongly deformed 
CNT-molccule heterostructureii^ Here, we present opti- 
mized EHT-parameters for bulk-silicon and benchmark 
its bandstructure against multiple target values such as 
bandedges and effective masses. We explore the trans- 
ferability of these EHT-parameters to different environ- 
ments by calculating the 2D-bandstructure for two re- 
construced silicon surfaces, silicon (100) (2x1) and (111) 
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(2x1), and compare them quantitatively to experiments 
and state of the art DFT-GW calculations. Finally, we 
use the silicon parameters to calculate the ID band dis- 
persion, density-of-states (DOS), and transmission (T) 
for a un-reconstructed silicon nanowire with and without 
hydrogen passivation demonstrating the capabilities of 
Extended Hiickcl Theory to model passivated surfaces of 
nanostructures. The EHT-paramctrization for silicon is 
shown to be quite transferable to different environments, 
so that a quantitative modeling of silicon-based devices 
becomes feasible. 

The paper is organized as follows: section II summa- 
rizes briefly the main features of Extended Hiickel The- 
ory. The optimized EHT parameters for bulk silicon 
along with the comparison of the bandedges and effec- 
tive masses to experimental target values are discussed 
in section III. We then investigate the transferability of 
the parameters by employing them to different silicon 
surfaces. Finally, the ID electronic structure of a silicon 
nanowire is determined including its surface passivation, 
and summarize our work in section IV. 



can be further generalized, so that heterogeneous sys- 
tems such as heterostructures and interfaces can be 
modclcdiiSiSSiSi The labels fj,, v refer to the atomic or- 
bitals, and is the overlap matrix between the orbital 
basis function (j)^ and <f> v , respectively. K e ht is an ad- 
ditional fit parameter usually set to 1.75 for molecules 
and 2.3 for solidsi22i2I Compared to the formal OTB- 
basis set, Slater basis functions are non-orthogonal, i.e. 
S ^ I, which provides an improved transferability of 
the model parameters with respect to changes in the 
environments^^^; the enhanced transferability can be 
justified by constructing orthogonal Lowdin orbitals from 
the non-orthogonal basis functions. Compared to the 
original AO's these Lowdin functions are more long- 
ranged to enforce orthogonality among different Lowdin 
functions over the entire domain. Consequently, the 
Lowdin basis becomes more sensitive to changes in the 
actual environment. 



II. APPROACH 



The silicon bandstructures for the bulk, the two re- 
construced surfaces, and the ID nanowire are calcu- 
lated within a non-orthogonal Slater-Koster scheme^ us- 
ing Extended Hiickel Theory to generate the overlap- 
and Hamiltonian-matrix elements S and H, respectively. 
Here, we briefly summarize the essential features of EHT, 
which is described in more detail in Refs>i242£. 

Extended Hiickel Theory is a semi-empirical theory to 
calculate the electronic structure of molecules and ele- 
mental solids. The most striking difference between EHT 
and orthogonal tight-binding (OTB) is that in EHT one 
works with explicit atomic-like orbital basis functions 
(AO), which are used to construct the matrix elements 
S and H. In turn, in orthogonal-tight binding the ba- 
sis functions are not known and used as a formal tool to 
construct all matrix elements of the Hamiltonian. The 
matrix elements are then usually adjusted to a reference 
bandstructure, for example. Compared to OTB in Ex- 
tended Hiickel Theory one adjusts only the diagonal ma- 
trix elements of the Hamiltonian (onsitc energies) and 
the parameters specifying the basis functions, which are 
Slater type functions (STO) i 18 ' 21 Since the basis func- 
tions are known, the overlap matrix S is calculated ex- 
plicitly and used to construct the off-diagonal matrix el- 
ements of the Hamiltonian (hopping) according toi^SS 




H^v — -^K e ht S^v {Hftft + H vv ) , (1) 

assuming that the Hamiltonian depends linearly on the 
overlap^ The original EHT-prcscription, cf. Eq. QJ, 



A concrete example where the transferability of pa- 
rameters becomes evident are structural deformations. 
For bulk-like systems empirical scaling rules have been 
developed^ so that effects of strain on the electronic 
structure can be studied using an orthogonal tight- 
binding approach employing the commonly used near- 
est neighbor approximation (NNA). Even though these 
empirical approaches for structural deformations are cali- 
brated and work well for bulk systems, the transferability 
of the bulk-optimized parameters along with the scaling 
rules have not been tested and benchmarked, for exam- 
ple, for reconstructed silicon surfaces. The difficulty in 
modeling (reconstructed) surfaces within OTB is that 
surface reconstruction 1). is usually accompanied by 
large structural changes beyond 2 — 5% and 2). is ac- 
companied by bond changes. With respect to point 1). 
it is not a priori obvious whether the scaling rules remain 
valid beyond small deformations for which they have been 
calibrated; regarding point 2). it is questionable if a 
nearest-neighbor approximation which works for bulk can 
be consistently applied to reconstructed surfaces. Using 
empirical scaling rules the hopping matrix elements be- 
tween neighboring atoms of the deformed structure are 
usually determined from the bulk hopping matrix ele- 
ment between the same atoms. In the case of the two 
dimer atoms the problem is that their initial bulk hop- 
ping matrix elements is zero, since their distance with 
f=a 3.8 Ais larger than the nearest-neighbor cut-off of 2.35 
A. However, to reproduce the experimentally observed pi- 
and 7r* surface bands it is important that the electronic 
structure model provides the correct hopping matrix el- 
ements between the two dimer atoms. 
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III. RESULTS FOR EHT-ELECTRONIC 
STRUCTURE FOR SILICON 

A. Si-Bulk 

To perform quantitative transport calculations 
through nanostructure materials, the free parameters 
of a semi-empirical tight-binding model have to be 
calibrated to experimental targets and/or bandstructure 
data obtained from other theoretical approaches. For 
Extended Huckel Theory Cerda and Soria have devel- 
oped EHT-parameters for several bulk crystal structures 
such as metals, semi-conductors, and compounds^ 
Specifically, for silicon these parameters have been 
optimized to match the bulk-dispersion of DFT-GW 
calculations of Rohlfing et al& at selected points within 
the 3D Brillouin zonei^i 

We start by optimizing the EHT-parameters using ex- 
perimentally determined bandstructure characteristics of 
bulk silicon such as bandedges and effective masses as 
targets^ We use the TBGreen codo^i to minimize the 
root-mean square error (RMS) between our EHT-bands 
and the targets via a conjugate gradient method as de- 
scribed in Refill. Since most of the targets refer to ex- 
periments done at low temperatures (5 — 10 KSi), we 
perform the minimization of the RMS-error at T = 0.0 
K. The 3D bandstructure of silicon is calculated using 
a sp 3 d 5 orbital basis for each silicon atom. In order to 
capture and to optimize the split-off gap Ao, the Hamil- 
tonian is spin-dependent through the spin-orbit coupling, 
LS P , used as additional fit parameter. 

Figure n shows the band dispersion for bulk-silicon us- 
ing the optimized EHT-parameters for a silicon atom^ 

In Table[I]we show the final results of our optimization 
along with their relative error (2nd and 3rd column) with 
the experimental target values as reference (3rd column). 
Our values for the bandedges as well as for the effective 
masses agree very well with the target values^ and are 
in their quality competitive to the state-of-the art or- 
thogonal sp 3 sd 5 s* tight-binding model for bulk-silicon^ 



B. Si-Surfaces for different Orientations 

We now investigate the transferability of our EHT- 
parameters for bulk-silicon and use them to calculate 
non-selfconsistcntly the surface bandstructure for two dif- 
ferent surface orientations. The semi-infinite surfaces are 
modeled by a finite slab consisting of a series of layers. 
Because of this truncation one introduces an additional 
free surface at the bottom of the slab and hence unphys- 
ical surface states, which can extend towards the true 
surface and thus tamper the surface band dispersion. To 
eliminate the dangling bond states at the bottom of the 
slab we passivate in each case the surface of (100) and 
(111) orientation by attaching hydrogen. 




l r x 

FIG. 1: Bandstructure of bulk-silicon calculated within EHT 
using the parameters given in Table— The Fermi level is at 
Ef = 0.0 eV as indicated by the horizontal solid line. The 
EHT parameters are optimized to experimental target values 
taken from Refi— ■ . 



To calculate the 2D silicon bandstructure of the recon- 
structed surface we use the unit cell coordinate of Ref^S 
for Si(100) (2x1) and for Si(lll) (2x1) from Ref^i as 
shown in Figurc[5] For Si (100) (2x1) we use a slab with 




FIG. 2: Structure of the unit cell for the two silicon 
surfaces!^!. Top:Si(100) (2x1) with the first 4 layers re- 
laxed and 9 bulk-like layers Bottom: Si(lll) (2x1) where the 
first 8 layers are relaxed and 12 bulk layers. In each case, 
the bottom of the surface is hydrogen passivated to remove 
dangling bond states. 

in total 13 layers, where the first 4 layers correspond to 
the reconstructed silicon surface, and the remaining 8 
layers correspond to positions of bulk silicon. The last 
layer is hydrogen to passivate the bottom of the slab. 
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TABLE I: Bandstructure characteristics for bulk silicon using Extended Hiickel Theory and spd-orbitals for each Si-atom. The 
EHT-parameters, cf. Table— have been optimized to experimental target values— (4th column). The last two column on the 
right are the fit values and errors based on an orthogonal tight-binding model using sp 3 d 5 s* orbitalsi- The effective masses at 
the L-valley are not well established, so that they are to strongly weighted in the optimization. 



The lattice vectors for Si(100) (2x1) are a x = 7.68A e x 
and &2 = 3.84A e y . Similarly, the Si(lll) (2x1) recon- 
structed surface contains 21 layers with the first 8 layers 
being relaxed, 12 bulk-like layers, and again the last layer 
hydrogen passivated. The 2D Bravais lattice vectors here 
are ai = 6. 65 A e x and &2 = 3.84A e y . 

Figure |3 shows the bandstructure of reconstructed sil- 
icon (100) (2x1) calculated within EHT (top) using the 
silicon parameters of Tabled. The band dispersion at 
the bottom of Figure |3| correspond to DFT-GW calcu- 
lations of Rohlfmg et alml with the it- and 7T*-bands in 
solid lines. As can be seen, our EHT-calculated n- and 
?T*-surface bands agree very well qualitatively in their 
shape with DFT-GW calculations as well as with PES 
experiments for the vT-band. 

Despite of the good qualitative agreement, which 
demonstrates the transferability of our parameters in 
capturing the essential physics of the surface bands, there 
are quantitative differences. The experimental indirect 
bandgap A ff *_ ff between the 7r* and 7r band, for example, 
is about s=s 0.8/0.9 eV, whereas the gap in EHT is with 
« 0.3 eV of similar order as in DFT-LDA calculations. In 
Table HT1 wc compare the bandgaps and bandedges calcu- 
lated in EHT at different points of the 2D Brillouin zone 
with those obtained by DFT-GW calculations^ and ex- 
periments. 

The differences between the theoretical approaches be- 
come more explicit in Figure [5] where only the E-k dis- 
persion for the 7T- and 7r*-surfacc bands is shown for the 
EHT- and the DFT-GW calculation of Rohlfmg et al. M 
The plots have been extracted from Figure |3] by digi- 
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TABLE II: Comparison of the bandgaps A and bandedges 
calculated in EHT, DFT-GWr^ and experiments for the sil- 
icon surface (100) 2x1 at different points of the 2D Brillouin 
zone. All values are in units of eV. 



tizing the respective tt- and 7r*-bands. The unoccupied 
7r*-band agrees quantitatively very well with the one ob- 
tained from DFT-GW over the entire Brillouin zone re- 
gion as shown in Figure^ Similarly, the 7r-band matches 
quantitatively DFT-GW as well over a wide range of the 
Brillouin zone, except within the T — J and the first 3rd 
of the J — K path. The latter domain includes the ir- 
maximum, which appears « 0.25/0.3 eV too high, so 
that the A w »_ ff bandgap is noticably underestimated. 
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FIG. 3: Surface bandstructure of unpassivated, reconstructed 
silicon (001) (2x1) calculated within EHT (top) using the 
parameters of Tabled The figure below is a DFT-GW 
calculation—. The two bands within the 2D-projected bulk- 
silicon bandgap correspond to the tt and n* bands of the sil- 
icon surface described by the ADM. Reprinted figure (mid- 
dle) with permission from M. Rohlfing, P. Kriiger, and J. 
Pollmann, PRB, 52, 1905 (1995). Copyright (1995) by the 
American Physical Society. 



The surface density of states (DOS) is shown in Fig- 
ure El The energy-resolved partial DOS is calculated 
for each dimer atom (upper and lower) and for the two 
deeper silicon layers away from the surface. The partial 
DOS of the upper dimer atom is located more closely to 
the valence band, whereas for the lower one it is near 
the conduction band, indicating that the 7t-surface band 
is formed from the upper-dimer atom, whereas the 7r*- 
band comes from the lower one^. Consistent with the 
7r-band dispersion of Figures |3] and 2] the PDOS of the 
upper dimer atom is too far away from the valence band, 
so that the ir* — it gap in the PDOS is too small. 

Away from the surface and approaching the bulk-like 
region the weight in the n- and 7r*-DOS decreases contin- 
uously (layer 4), and completely disappears once a deeper 
bulk- like layer is reached (layer 8), so that the original 
bulk-bandgap of w 1.2 eV is recovered. 




FIG. 4: Comparison of the n- and 7r*-surface bands calculated 
within EHT and DFT-GVv^. The dispersion of the tt- and 
7r*-bands of Figure |H] have been digitized. The data of the 
red curve are adapted with permission from Figure 3 of M. 
Rohlfing, P. Kriiger, and J. Pollmann, PRB, 52, 1905 (1995). 
Copyright (1995) by the American Physical Society. 
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FIG. 5: Density-of-states of un-passivated silicon (100) (2x1) 
surface calculated for different layers starting from the bulk- 
like 8th layer (top) towards the (100) surface consisting of the 
two dimer atoms. The two peaks in the DOS at each dimer 
atom correspond to the n- and n* bands. 
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As second example, we look at the surface bandstruc- 
ture of reconstructed Si(lll) 2x1 as shown in Figure 
using the EHT-parameters in Tabled. Similar to the 
previous case, the overall shape of the tt- and 7r*-surfacc 
bands match qualitatively well with DFT-GW calcula- 
tions of Rohlfing et alfiK In Tabic Ell we compare our 
EHT-bandedges and gaps (1st column) at two specific 
points J and K of the 2D-Brillouin zone with DFT-GW 
calculationo 35 i 36 and PES/IPES experiments. The values 
for the bandedges as well as for the gaps agree quantita- 
tively well among all three calculations, and show also a 
good agreement with PES/IPES experiments, where the 
error in the energy resolution is typically 150 — 200 meV 
depending on temperature and incident energy of the 
electrons. A more extended comparison with PES/IPES- 
cxperiments turns out to be very limited, since the tt- and 
the 7r*-bands are not as well experimentally determined 
as in the case of silicon (100) (2x1). Contrary to the pre- 



entire range of the Brillouin zone as shown in Figure [7| 
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FIG. 6: tt- and tt* surface bandstructure of the unpassivated, 
reconstructed silicon (111) (2x1) calculated in EHT (top). 
The figure at the bottom shows the dispersion calculated 
within DFT-GGA of Rohlfing and Louie^, respectively. The 
bottom figure is reprinted with permission from M. Rohlfing 
and S.G. Louie, Phys.Stat. Solidi (a), 175, 17 (1999). Copy- 
right (1999) by the American Physical Society. 

vious case of Si(100) (2x1), we find for Si(lll) (2x1) that 
both tt*- and 7r-band calculated in EHT agree quantita- 
tively very well with DFT-GW calculations of Northrup 
et td~- and in particular with Rohlfing et al~ over the 
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FIG. 7: Comparison of the tt- and 7r*-surface bands calculated 
within EHT and DFT-GW^^. The dispersion of the tt- and 
7r*-bands of Figure |3] have been digitized. The data in the 
green curve are adapted with permission from Figure 2 of J.E. 
Northrup, M.S. Hybertsen, and S.G. Louie, Phys. Rev. Lett., 
66, 500 (1991). Copyright (1999) by the American Physical 
Society. The data in the red curve are adapted with permis- 
sion from Figure 2 of M. Rohlfing and S.G. Louie, Phys.Stat. 
Solidi (a), 175, 17 (1999). Copyright (1999) by the"American 
Physical Society. 





EHT (Tabled) 


DFT-GW^ 


DFT-GW^ 


Exp. 


Aj 


0.66 


0.60 


0.68 




A K 


0.91 


0.82 


0.92 




Kj 


-0.05 


-0.07 


0.0 


0.09 


Ef 


0.61 


0.53 


0.69 


0.67 


E\ 


-0.11 


-0.04 


-0.04 




Ek 


0.80 


0.78 


0.88 





TABLE III: Comparison of the bandgaps A and bandedges 
(in units of eV) calculated in EHT, DFT-GW2^, and exper- 
iments for the silicon surface (111) (2x1) at different points 
of the 2D Brillouin zone. 

In the two previous cases we explored the transfer- 
ability of the EHT-parameters, cf. Tabled , optimized 
for bulk-silicon by applying them to other environments 
such as re-constructed surfaces for silicon (100) (2x1) 
and (111) (2x1). Without any re-parametrization the 
experimentally observed tt- and 7r*-surface bands are re- 
produced qualitatively in their overall shape, and in the 
case of silicon (111) (2x1) we also achieve a good quan- 
titative match compared to PES/IPES experiments and 
DFT-GW calculations as well. As discussed, quantitative 
differences exist, particularly for the indirect bandgap 
Att.-tt for silicon (100) (2x1), which is underestimated 
similar to DFT-LDA pseudopotential calculations^ 

One reason for the obvious quantitative differences, 
particularly the wrong position of the 7r-band above the 
valence might be due to the non-self consistent calcula- 
tion of the bandstructurcs for the reconstructed silicon 
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(100) and (111) surfaces. Calculating the total non-self 
consistent charge on each dimer atom by integrating the 
LDOS gives a charge of 4.13e on the upper and 3.75e on 
the lower dimer atom. The total charge of the two dimer 
atoms is about 7.88e. In turn, a self-consistent calcu- 
lation of the dimer atom charge using SIESTA(2L2& re- 
sults in a total charge on the upper dimer atom of about 
4.0c, whereas the lower one has 3.88e. Note, that in 
the EHT non-scf as well as in the SIESTA scf-case the 
total charge on the two-dimer system is the same and 
is effectively positive with respect to their total num- 
ber of valence electrons of 8. Qualitatively, what one 
would expect is that the upper dimer atom which car- 
ries initially too much negative charge, looses parts of 
it under self-consistency and charge is transfered partly 
to the lower dimer atom making, which becomes in turn 
more negative after self-consistency. The overall effect 
of self-cinsistency is then mainly to redistribute charge 
among the two dimer atoms; the initial 7r-band which 
consist of the upper (less negative) dimer would float 
down, whereas the 7r*-band would float up since the lower 
dimer atom is more negative. Both bands float in op- 
posite directions, so that the indirect bandgap A^*-^ 
increases. 



C. Si-Nanowire (100): H-Passivation in EHT 

As final example, we use the silicon parameters, cf. 
Tabled in combination with the ones for hydrogen to 
demonstrate how to passivate surfaces of nanostructures 
in EHT by means of a silicon nanowirc. We use a wire 
along the (100) direction with rectangular cross section 
and sidelength D — 1.5 nm as shown in Figure |H1 For 
simplicity we assume that the wire is un-relaxed, i.e. the 
positions of the silicon atoms of the wire correspond to 
the positions of bulk silicon. For the surface we consider 
two cases as illustrated in Figure |H1 i) a wire without H- 
passivation and ii) a wire where we have explicitly added 
hydrogen atoms to saturate the dangling bonds. 

In Figure[31the ID-dispersion relation is shown for the 
two structures, cf. Fig. |H1 As can be seen in the left 
part, the "bulk" -like bandgap of the silicon nanowire is 
covered by dangling bond states due to the unsaturated 
bonds of the surface atoms. The entire bands of the 
dangling bonds are completely removed after the wire is 
passivated by explicitly attaching hydrogen (bottom), so 
that the bulk-bandgap of w 2.1 cV is recovered. Note, 
that for the case of a H-passivated wire we also assume 
that the atoms remain in their bulk-positions. 

The approach in Extended HuckcL Theory to passi- 
vate the SiNW surface (and possible other surfaces) by 
physically attaching hydrogen atoms, for example, differs 
from the methods described in Ref^ within an orthogo- 
nal tight-binding scheme. In the latter case, the dangling 
bonds are removed by transforming first the Hamiltonian 
from a \s), \p x ), \Py), \Pz) representation to a represen- 
tation using sp 3 -hybrid basis functions. One then has to 




FIG. 8: Front (top) and sideview (bottom) of a silicon 
nanowire with rectangular cross section and sidelength D = 
1.5 nm along the (100) direction. For the un-passivated wire 
(left) the unit cell contains 81 silicon atoms, and in the H- 
passivated case (right) the total number of atoms is 117. 




FIG. 9: ID bandstructure and transmission per spin of an 
un-relaxed silicon nanowire along the (100) wire axis. The 
E — k and transmission on the left are for the unpassivated 
surface, and on the right the surface is hydrogen passivated to 
remove the dangling bond states. The transmission shown in 
the insets takes integer values where each channel contributes 
one unit quantum conductance Go = e 2 /h per spin. 



check which bonds of the surface atoms are not saturated 
and to increase then the respective orbital energy of the 
dangling bond. In raising the orbital energy manually 
one formally mimics the effect of the hydrogen atoms 
onto the silicon surface atoms as if hydrogen were at- 
tached. The amount, however, by which the dangling 
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bond state energy has to be shifted is a priori not known 
and has to be determined empirically^. 

In turn, in Extended Hiickel Theory the surface passi- 
vation is controlled by the passivation atoms, which are 
an explicit part of the entire structure. Once the struc- 
ture is specified, the amount by which the silicon levels 
are shifted is determined by the chemical species of the 
passivation atoms. The hybridization causing the shift 
of the silicon dangling bond states is naturally incorpo- 
rated within the Hamiltonian matrix through the EHT- 
prescription, cf. Eq.Q. 

IV. SUMMARY 

We applied Extended Hiickel Theory (EHT) to de- 
scribe the electronic structure of several silicon-based 
structures. A quantitative benchmark of the bulk-silicon 
electronic structure was achieved by optimizing the EHT- 
parameters to experimentally determined target values 
specific for the bulk-bandstructure of silicon. Our de- 
veloped parameter set for silicon provides a quantitative 
good agreement with these targets and arc competitive 
with the results based on orthogonal-tight binding for 
silicon. 

The transferability of the parameter set was investi- 
gated by applying them to different environments such 
as reconstructed silicon for two different surface ori- 
entations. Within our EHT-approach using the bulk- 
optmized parameters we could reproduce the essential 
features of the surface band dispersion, in particular the 
7r- and 7r*-surfacc bands experimentally well established. 
A quantitative comparison of our EHT-bandstructure 
with DFT-GW calculations as well as PES/IPES experi- 
ments shows qualitative and quantitative agreement, par- 
ticularly for silicon (111) (2x1) surface. However, our in- 
direct bandgap is lower compared to PES/IPES 
experiments and of the same order as DFT-LDA calcu- 
lated gaps. This discrepancy is partly due to the non-self 
consistent calculation of the surface electronic structure 
of reconstructed silicon, where the charge redistribution 
due to structural changes is discarded. We expect that a 
full 3D self-consistent solution of the electronic structure, 
for example within a Complete iVeglcct of Differential 
Overlap scheme (CNDO), can correct for the discrepan- 
cies, which would further increase the transferability of 
our EHT-parameters for silicon. 

Using the silicon nanowire as example we used the 
EHT-parameters for silicon and demonstrated a generic 
approach to surface passivation of nanostructure materi- 
als by physically attaching hydrogen atoms to the SiNW 
surface. The respective dangling bond states are removed 



from the bandgap region of the wire in a systematic man- 
ner without the need to shift dangling bond molecular 
levels by hand. 

The flexibility of EHT demonstrated here and in our 
previous paper (Part I)i& opens the door to study elec- 
tronic structure and transport through molecular het- 
erostructures as well as larger nanostructures preserving 
the atomistic features of the system. We believe, that 
Extended Hiickel Theory is a good practical comprimisc 
between rigorous, but computationally expensive DFT- 
based approaches, and orthogonal tight-binding meth- 
ods, which might not be suitable for large structural de- 
formations beyond 2 — 3%. More coarse-grained models 
such as effective mass might be even prohibitive for the 
same purpose, since they inherently fail to account for 
bonding. The main appeal of EHT is that it does capture 
bulk- as well as surface physics along with bonding chem- 
istry at heterointerfaccs and molecular hctcrostructures - 
including large structural deformations - all within a uni- 
fied semi-empirical frameworki&SiS. Silicon nanowires, 
for example, are interesting candidates for new channel 
materials for novel MOSFETs. A widely discussed and 
still open question is whether structural relaxation of the 
wire significantly affects its electronic structure as well as 
the overall device performance. Intuitively, one expects 
that the electronic properties particularly of small diam- 
eter wires are dictated by their surface. Some aspects 
of this problem is currently investigated where we em- 
ploy Extended Hiickel Theory to describe the electronic 
structure of the relaxed wire for different diameters and 
orientations4i Extended Hiickel Theory as it stands also 
needs further improvement to fully utilize its capabili- 
ties and to establish it as a methodological tool towards 
a quantitative modeling of quantum transport through 
nanostructures. This is left for future work. 
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